Vortex Lattices in Stripe Domains of Ferromagnet/Superconductor Bilayer 
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The continuum theory of domain structures in ferromagnetic/superconducting bilayers fails when 
the equilibrium domain size becomes comparable with effective penetration depth A. Instead, a lat- 
tice of discrete vortices must be considered. Here, we report our results on the discrete vortex lattices 
in stripe domain structures of ferromagnetic/superconducting bilayers. The vortices are assumed 
to be situated periodically on chains in stripe domains. We study the configurations containing up 
to three chains per domain, and calculate their equilibrium energies, equilibrium domain size and 
vortex positions through a method based on London-Maxwell equations. In equilibrium, the vortices 
in the neighbor domains are half-way shifted while they are next to each other in the same domain. 
Additionally, more vortex chains per domain appear spontaneously depending on magnetization and 
domain wall energy. 



PACS Number(s): 74.25.Dw, 74.25.Ha, 74.25. Qt, 74.78.-w 



I. INTRODUCTION 

Heterostructures made of superconducting (SC) and 
ferromagnetic (FM) pieces not only give rich physical 
effects that are not observed in individual subsystems, 
but also offer new devices that can be tuned by weak 
magnetic fields [1,2]. One of the realizations of such 
heterostructures is ferromagnetic/suprconducting bilayer 
(FSB). Earlier Lyuksyutov and Pokrovsky noticed [3,4] 
that in a bilayer consisting of homogeneous SC and FM 
films with the magnetization normal to the plane, SC 
vortices occur spontaneously in the ground state, even 
though the magnet does not generate a magnetic field in 
the SC film. 

In previous work, we presented a theory of such vortex- 
generation instability and the resulting vortex structures 
[5] . We showed that due to this instability domains with 
alternating magnetization and vortex directions occur in 
FSB. In that study, we treated these domain structures in 
the continuum regime in which the domain size is much 
larger than the effective penetration depth, which is de- 
fined to be A = X 2 /d sc , where the London penetration 
depth A is much larger than the thickness of supercon- 
ducting film d sc [6]. Under the continuum aproximation, 
we found that the energy of stripe phase was minimal. 
The equilibrium domain size and the equilibrium energy 
for the stripe structure was found as [5] 
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where m = m — £ v /<f>o, £v = (0o/167r 2 A) ln(A/£) is self- 
energy of a vortex, Edw is domain wall energy per domain 
wall length, .4 is the domain's area and C ~ 0.577 is the 
Euler-Mascheroni constant. If £dw < 4m 2 , the contin- 
uum approximation becomes invalid, since L eq becomes 
on the order of or less than A (see Eq.(l)). However, 
it can be recovered by considering the discrete lattice of 
vortices instead. 

In this paper, we study the discrete latice of vortices in 
stripe domain structures in FSB via a method based on 
London-Maxwell equations, which is developed elsewhere 
[7] . The extension of the method to periodic systems is 
also introduced for the case of square magnetic dot arrays 
on a SC film [8]. Here we adapt it for the discrete vor- 
tex lattices in SC/FM bilayers. In doing so, we assume 
that vortices and antivortices sit periodically on chains in 
the alternating domains of magnetization and vorticity. 
The problems we consider here are; i) how the vortices 
and the antivortices are positioned on the chains; ii) how 
the equilibrium domain size changes, depending on the 
magnetization and the magnetic domain wall energy in 
the presence of the vortices; iii) if the spontaneous ap- 
pearence of domains structures with different number of 
vortex chains is possible. In order to solve these prob- 
lems, we first propose five different configurations of the 
vortex and the antivortex chains, in which at most two 
chains per stripe is considered. Next, we calculate their 
equilibrium energies by means of numerical methods and 
find the most favorable case among them. Our calcula- 
tions show that in equilibrium structure, vortex chains 
are half shifted in the adjacent domain while they are 
next to each other in the same domains. Inspired from 
this result, similar configuration for three vortex chains 
per domain is considered, and its equilibrium energy is 
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calculated for various values of magnetization and do- 
main wall energy. Comparison of equilibrium energies of 
cases with one, two and three chains per domains shows 
that at lower values of magnetization and domain wall 
energy case with two chains is favorable whereas , three 
chain cases wins at higher values of magnetization and 
domain wall energy. Additionally, single chain case does 
not win over the ones with two and three chains per do- 
main under any circumstances. 



The outline of this paper is as follows: In the follow- 
ing section, we present the method for discrete case and 
its application to configurations with a single and double 
vortex chains per domain. In the third section, we briefly 
present our results on the proposed five configurations. 
The fourth section is devoted to the study of the case 
with three chains per domain. The last section consists 
of the conclusions and discussion. In the appendix, we 
give the details of the methods and mathematical tricks 
in series calculations. 



(a) 1st Case 



(b) 2nd Case 
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II. METHOD 



(c) 3rd Case:l 



(d) 3rd Case:2 



In this section, we introduce a method based on the 
treatment of vortices in the discrete lattice. We study the 
lattices of discrete vortices only in the stripe phase. In 
the continuum approximation, it is found that the vortex 
density increases at the closer distances to the magnetic 
domain walls. Based on this fact and the symmetry of the 
stripe domain structure, it is reasonable to consider that 
the vortices and antivortices form periodic structures on 
straight chains along the y direction. Even though it is 
not clear how many chains are associated with each do- 
main, we can still make progress toward understanding 
discrete vortex lattices. To this end, five stripe domain 
configurations in which vortices are situated periodically 
on chains are proposed. From this point on, the con- 
figurations with N vortex chains per stripe domain are 
labeled as N state. 

In two of proposed cases, there is one chain per stripe 
(N = 1 states), located in the middle of the domain. In 
this case, two configurations of vortex lattice are possi- 
ble; first, the vortices and the antivortices in a neighbor- 
ing domains are alongside to one another (see Fig. 1(a)), 
second, they are shifted by half period 6/2 along the y 
direction, where 6 is the distance between two nearest 
vortices on the chain (see Fig. 1(b)). 



+ + 



+ + 



+ + 



(e) 4th Case 



(f) 5th Case 



FIG. 1. Proposed configurations for N 
states. 



1 and N 



In the other three candidate lattice structures, there 
are two chains per stripe domain N — 2 states, at the 
distance a from the magnetic domain walls. The possi- 
ble situations in the two chain cases are as follows; First, 
chains in the same domains are shifted by a half period 
6/2 along y direction (see Fig. 1(c)) . However, the neigh- 
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bor vortices and antivortices are next to each other. Sec- 
ond, they are just shifted by a half period 6/2 (see Fig. 
1(e)). Third, the vortices and the antivortices are simply 
side by side on the chains (see Fig. 1(f)). 

Our next step is to write the energies of these five can- 
didates. To this end, we use the energy equations for 
periodic systems developed elsewhere [8]. 
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In those equations, the vortex configurations differ by 
their form- factors. We can obtain them from Fa = 
J2r nie lG ' ri , where the G's are the reciprocal vectors 
of the periodic structures, the are the positions of the 
vortex centers, and Hi are the charge of the vortex. In 
our proposed models, G = ((2r+ 1)£, 2s£) and rn = ±1. 
Table I gives the form factors of each configuration in the 
order they are described above. 

TABLE I. The form-factors of vortices in the proposed con- 
figurations. 
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Note that, in Table I, the form-factor for the third con- 
figuration also belongs to the case in which the vortex 
and the antivortex chains are shifted by half period only 
in the neighbour domains, not in the same domain (see 
Fig. 1(d)). Since information about the vortex lattice 
is carried only by the form factors, there is no need to 
consider the above-mentioned case separately. 

In our calculations, the divergent part of the series 
must be extracted carefully. We show below the detailed 
analysis of series equations for each candidate. We start 
with the self interaction energy of the magnetic layer 
Umm, since it is the same for each configuration. For 
the periodic structures, it is given by (5). Direct sub- 
stitution of the Fourier coefficient of the stripe phase 
m Z G = ff (2r+i) mto Eq.(5) gives the self-interaction of 
the magnetic layer per unit cell as 
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where (x) is the polygamma function of zeroth order 
[9]. In our numerical calculations, we write the logarith- 
mic term in (6) as ln(A//)+ln(L/A) and then incorporate 
the —4m 2 ln(A/Z) term in the renormalized £ r J^ ■ Another 
energy term with a divergent series is the vortex energy, 
in general given by (3). The logarithmic divergence in 
this term stems from the vortex self-energies. We first 
split (3) into two parts as follows: 
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Note that the area of the unit cell is 2Lb. The first term 
of the series above contributes to the self-energies of the 
vortices; whereas, the second term is the vortex-vortex 
energy and will be left in the series form. For each form- 
factor in Table I the series in the first term can be trans- 
formed to the form of 

Y^T=-oo X^-oo l/((2r + l) 2 x 2 + s 2 ) where x is constant, 
and depends on the form-factor. A detailed analysis of 
such series is given in Appendix A. 

The next step is to find the vortex energy and the 
interaction energy of the magnetization and vortices 
for each configuration. In the calculation of 

u-mv j we 

take the Fourier coefficient of the magnetization to be 
(2r+i) $(Gy)- The fact that the stripe is infinite along 
the y direction results in the additional term 2ir5(G y ). 
However, it does not play any role in the calculation of 
u mm - For numerical analysis, these energies must be ex- 
pressed in terms of dimensionless parameters. To this 
end, we define dimensionless variables A = A/ L , b = bj L 
and idw — e d«T^/ e o- The total energy U is measured in 
units of Eq/ A 2 . In addition, we introduce the dimension- 
less magnetic energy as U m m — «mm/(£o/A 2 ). In terms 
of these parameters, the energy of the first configuration 
reads 



3 



um -tMM + 2tf ,<A) - ^ - ^ <A) ) + °- + 



where, 



~ coth((2r +!)%)-! 



/ (1) = V 

J VV / j 



V.S — — 00 



((2r + l) 2 + ^)(1 + 2^(2^+1)2 + ^)) ' 



" ^ (2r + l)(l + 27rA(2r + l))" 1 ; 



The form factor for the second configuration survives of the second configuration is found to be 
only for even values of s. Then, the dimcnsionlcss energy 

26 V \e G A£/ &7T £o / 



where /i 2 ] = f£l and, 



(2) ~ coth((2r + l)^)-l 
2r + 1 

OO 

f(2) = 

J VV 



% -V ± ( 1 1 ) 

— oo ((2r + l) 2 + i^)(l + 2^^+1)2 + iff)) ' 



In the third configuration, as in the second configu- 
ration, only even values of s contribute to the energy. 
In the first two configurations, the square of their form- 
factors enter the vortex energy as a constant. However, 
in this case, the square of the sine function appears. In 



Appendix A, it is shown how to calculate the series in 
the presence of such functions. Introducing the dimen- 
sionless parameter a = a/L, the energy functional of the 
third configuration becomes 
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In the fourth configuration, the square of the form- of s give different contributions. Then, we can calcu- 

factor is: late the vortex energy for even s and odd s separately. 

|Fg| 2 = 2 — 2(— l) s cos((2r + l)7ra). Even and odd values Employing similar techniques, we find 

^ (4) = € ( ln (-On) + 2 ^ 4) ( A < 5 ) - r ^ )(A ' a ' *) - — /^( A > «)) + #™» + ^ ( 14 ) 
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The form-factor for the fifth case resembles that of the integers contributes to the sum over s in the fifth case, 

third case with an exception. That is, in the third case, Keeping this in mind, it is straightforward to obtain the 
only even values of s are taken into account, while all dimcnsionless energy for the last case as 

t/ (5) = | ( J^) - ln(cot ( ™) + 4/<»>(A,a) - £/<?(A,«,6) - ^/£>(A,a) 

+ U mm + £dwK (16) 
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III. RESULTS FOR N = 1 AND N = 2 STATES 

In this section, we present our results based on nu- 

(i) 

merical calculations. The series in fi, converges very 
fast when r max > 200 , while the series in f^v and con- 
verges rather slowly, although, when r max > 4000 and 
Smax > 4000, the results do not change up to 6th deci- 
mal point in the energy, where i labels the particular do- 
main configuration. To make sure of this accuracy in the 
calculations, we take r max — 600 for fy l \ r max = 5000 
for f$ v , and r max = 5000 and s max = 5000 for f$ in 
Eqs.(8, 10, 12, 14, 16). In addition, the respective error 
deviations for X/L,b/L and a/L in numerical calculations 
are ±0.005,±0.0005 and ±0.0125. 

In the numerical minimization of Eqs.(8, 10, 12, 14, 
16), we take ln(4A/(e c £)) = 5.57. Changing m0 o /eo 
at fixed idw, we calculate the minimal energy of each 



configuration. We first investigate when these configura- 
tions become energetically favorable in the system. To 
this end, we check where the equilibrium energies of the 
configurations first become negative. 

In our analysis, we first identify two regimes of inter- 
est: discrete regime ( Edw < 4m 2 ) and continuum regime 
( Edw > 4m 2 ). Using our parameters idw and m<po/eo, 
the inequality for discrete regime can be expressed as 

£dw - — \-^Hf)- 1 ) ' (18) 

From Eq.18, the minimum value of m4>o/so for different 
values of idw can be determined. These values are given 
in Table.II 

TABLE II. idw versus {m<f>o/£o) m in- The column on the 
left is input. 
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(m0o/eo)min 


0.01 


6.198 


0.10 


7.557 


1.00 


11.853 


10.00 


25.439 



When m4>o/so is greater than its minimum value for 
a certain id w , the system is then in the discrete regime. 
Next, we analyze the proposed cases for N = 1 and N ~ 2 
states according to numerical values given in Table. II. 
The equilibrium of energies for these cases in discrete 
and continuum regimes are given in Table. III. 



TABLE III. Equilibrum energies for proposed configurations. Two columns on the left are input. 



£ dw 


(rn(f>o/Sv)min 


Ui 


Ul 


Us 


Ua 


Us 


0.01 


5 


-2.58454176 


-2.58455668 


-3.36407195 


-3.36404625 


-3.36404615 


0.01 


20 


-65.98296440 


-65.98296500 


-89.20105311 


-89.20030961 


-89.20030943 


0.1 


5 


-2.54211623 


-2.54211637 


-3.33057991 


-3.33054828 


-3.32949063 


0.1 


20 


-65.87136440 


-65.87136500 


-89.10475311 


-89.10310961 


-89.10310943 


i 


5 


-2.16635495 


-2.16637343 


-3.01972350 


-3.01972355 


-3.01972347 


i 


20 


-64.78187766 


-64.781878054 


-88.15826391 


-88.15826394 


-88.15826367 


10 


5 


-0.35644861 


-0.35667498 


-1.21258278 


-1.21232397 


-1.21232396 


10 


25 


-95.00241100 


-95.00247780 


-134.28090345 


-134.27589780 


-134.27589769 



G 



In our numerical calculations, we find that all pro- 
posed configurations are stable both in discrete and con- 
tinuum regimes. This indicates that our method works 
well in both regimes. Our numerical calculations also 
show that the third configuration wins over the other 
cases. Nonentheless, this is not enough information for 
us to understand the equilibrium structure, since 

[/(3) 

corresponds to two different cases with the same struc- 
ture factor. At this point, we need further analysis to 
determine which configuration is more likely. This can 
be done from simple physical considerations. Namely, 
in FSB, the equilibrium structure is determined by the 
competition between vortex-vortex interaction and vor- 
tex - magnetization interaction. The former favors vor- 
tices and antivortices in neighbor domains to line up in 
transverse direction ( perpendicular to magnetic domain 
wall ), whereas the latter prefers vortices and antivor- 
tices to be shifted so that gain in energy is maximized. 
When vortices are next to each other on the either side of 
the magnetic domain wall, the magnetic fields they pro- 
duce cancel out each other. From the numerical results, 
it is obvious that vortex-magnetization interaction wins 
the competion, and results in half-way shifting of vor- 
tices, if one compares the energies of 1st and 2nd cases. 



Then, vortex-magnetization interaction is the dominant 
factor. By the same token , one can understand what 
is going on in double vortex chain configurations. For 
instance, in the fifth configuration, energy gain due to 
vortex-magnetization interaction is diminished, since all 
the vortices are side by side. This explains why equilib- 
rium energy of the fifth configuration is higher than those 
of third and fourth cases. In the fourth configuration, 
the vortices and antivortices in the neighbor domains 
are half-way shifted, so that this configuration must be 
preferred over the one in which they sit side by side in 
the neighbor domains according to the above arguments. 
However, alternative configuration for third case has two 
chains half-way shifted in the neighbor domains instead 
of one chain as in the fourth case. Therefore, one might 
expect gain is even more than that in the fourth con- 
figuration. Another interesting result is that the system 
does not favor N = 1 state at all. Actually, this docs 
not surprise us since, in the continuum approximation, 
we found that the vortex density increases near the mag- 
netic domain walls. This fact already suggests that the 
system favors vortex chains being near the magnetic do- 
main walls rather than a single chain in the middle of the 
domain. 
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(a) 1st Casc:L/A (b) 2nd Case:L/A 




(c) 3rd Case:L/A (d) 4th Case:L/A 




o 1 1 1 1 1 1 

5 10 15 20 25 



/e 

(e) 5th Case:L/A 



FIG. 2. L/A versus m0o/eo for = 1 and N = 2 states. 
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e dw = 001 
E d „ = 0.1 



10 15 20 



10 15 



Edw - 0.01 
E dw = 0-1 

£ dw = 1 
EdW = 10 



(a) 1st Case:b/A 



(b) 2nd Case:b/A 




E dW = 001 
EdW " 0.1 
E d w= 1 
E d w= '0 




(c) 3rd Case:fe/A 



(d) 4th Case:b/A 




FIG. 3. b/A versus m^o/eo for N = 1 and N = 2 states. 




FIG. 4. a/A versus m0o/eo for N = 1 and N = 2 states. 



In numerical calculations equilibrium domain size L/A, 
vortex- vortex distances on the same chain b/A and 
vortex-magnetic domain wall distances a/ A are also cal- 
culated. These results are depicted in Figs. 2, 3, 4. At 
fixed Edwi the further increase of m<po/ea shrinks the 
domain width, while higher domain wall energy favors 
larger domain width at fixed m0o/£o, as in usual ferro- 
magnets. That is to say, ferromagnet favors narrower do- 
mains to minimize the demagnetization energy whereas 
domain wall energy makes domains wider. The compe- 
tition between these two energies determine the domain 
size. Here, the parameter m^o/eo plays the role of de- 
magnetization energy. Domain wall energy does not af- 
fect the distance between the vortices located on the same 
chain. However, at larger values of m</>o/£o the vortices 
on the same chain gets closer. This implies that the unit 
cell area shrinks , and consequently vortex density per 
area increases. Results of unit cell areas and vortex den- 
sities are discussed with details in the next section. 



IV. N = 3 STATE 



In this section, our aim is to understand whether the 
domain structures with different number vortex chains 
appear spontaneously or not. To answer this question, 
we consider only N = 3 state and compare its equilib- 
rium energy with that of N — 2 state for various values 
of magnetic domain wall energy and magnetization. We 
first show the derivation of energy equation for N = 3 
state. Following the physical arguments in the previous 
section, we consider that the equilibrium configuration 
is superposition of the 2nd case and alternative 3rd case 
(see Fig. 5). 
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FIG. 5. The vortex lattice in the first configuration. 



Then, its structure factor equals to sum of structure 
factors of those cases, which is given by 



F G = (2ism(G x a) + i(-l) r )(l + (-l) s ). 



(19) 



As seen in second section, self-vortex energy and vortex- 
vortex interaction depend on |-Fg| 2 , whereas Fg is used 
in vortex-magnetization interaction energy. Then, it is 
easy to see vortex-magnetization energy in this 
the sum of those of 2nd and 3rd case. We then calcu- 
late only the vortex energy here. The square of struc- 
ture factor is I^gI 2 = l^l 2 + \F^\ 2 + {F^)*(F^) + 
(Fq^)*(F(?). The first two terms give the vortex energy 
contributions of 2nd and 3rd case. Then, it is necessary 
to calculate contribution from the cross term only, which 
reads 



4A 2 ^ /cot((2r + 1W6/4) - l\ 2A 2 , lA ... 

U C ross = -g-5Zsin((2r + l)7ro)(-l) r ( U 2r l 1 J +— ln|cot(7r/4-7ra/2)| 



(20) 



2A 2 



2 b 2 



E 



sin((2r + l)7r5)(-l) r 



(2r + l) 2 + 16s 2 /6 2 i + 27 rA v /(2r + l) 2 + 16s 2 /6 2 ' 



Then, total energy becomes 

r(2) 



( <"= 3 - i (^) + ^ - ™ - ^< A >) < 21 > 

Ai ^ - In(cot(7ro)) + 4/< 3 )(A, 5) - £/$(A, 5, 6) - 1^/( 3 )(A, a)^ + U cross + U mm + e dw A. 



Note that U mm is described in section II. Next, Eq.22 
is minimized w.r.t. A, b and a for different values of idw 
and mcf>o/so. In numerical calculations, the series are cal- 
culated as described in the previous section. The com- 
porison of equilibrium energies of ./V = 1, N = 2 and 



N = 3 states at fixed values of e dw and m<j)Q /sq are given 
in Table. IV. Note that N = 1 state represents the second 
configuration in section II, since it is the most energet- 
ically favorable configuration between two single vortex 
chain configurations. 



TABLE IV. Equilibrium energies of two chain and three chain states. Vortex density of unit cell n equals N ce u / S e i , where 
N ce ii = 27V is the number of vortices in a unit cell and S e i = 2Lb is equilibrium unit cell area. The two columns on the left are 
input. 
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TTld)n / £n 


jJN=l 




n 


jJN=2 




n 


V N=3 




71 


0.01 


5 


-2.5846 


2.70 


0.74 


-3.3641 


4.11 


0.97 


-3.2576 


6.64 


0.90 


0.01 


10 


-13.9742 


0.73 


2.74 


-18.5183 


1.04 


3.82 


-18.8330 


1.58 


3.79 


0.01 


15 


-35.0799 


0.36 


5.56 


-47.0362 


0.50 


8.02 


-48.7311 


0.75 


8.02 


0.01 


20 


-65.9830 


0.22 


8.98 


-89.2011 


0.30 


13.24 


-93.4879 


0.44 


13.49 


0.1 


5 


-2.5421 


2.75 


0.73 


-3.3306 


4.22 


0.95 


-3.2324 


6.64 


0.90 


0.1 


10 


-13.9015 


0.74 


2.71 


-18.4598 


1.08 


3.71 


-18.7780 


1.58 


3.79 


0.1 


15 


-34.9848 


0.36 


5.52 


-46.9570 


0.50 


8.02 


-48.6699 


0.75 


8.02 


0.1 


20 


-65.8714 


0.23 


8.74 


-89.1048 


0.30 


13.24 


-93.4122 


0.44 


13.49 


1 


5 


-2.1664 


3.39 


0.59 


-3.0197 


4.84 


0.83 


-2.9941 


7.13 


0.84 


1 


10 


-13.2169 


0.82 


2.44 


-17.8987 


1.13 


3.54 


-18.3421 


1.68 


3.56 


1 


15 


-34.0699 


0.39 


5.14 


-46.1860 


0.52 


7.67 


48.0673 


0.77 


7.78 


1 


20 


-67.7819 


0.24 


8.40 


-88.1583 


0.31 


12.79 


-92.6655 


0.45 


13.18 


10 


5 


-0.3567 


14.2 


0.14 


-1.2126 


10.99 


0.36 


-1.4638 


14.16 


0.42 


10 


10 


-8.6308 


1.58 


1.27 


-13.6758 


1.78 


2.24 


-14.8269 


2.37 


2.53 


10 


15 


-27.1877 


0.63 


3.18 


-39.9099 


0.74 


5.41 


-42.8834 


0.99 


6.05 


10 


20 


-56.0242 


0.34 


5.80 


-80.1062 


0.41 


9.66 


-86.0271 


0.56 


10.78 



As seen from Table. IV, at low magnetization and small 
domain wall energies, N — 2 state is energetically favor- 
able. However, N = 3 state wins over at high values 
of m^o/^o and idw This picture resembles the one we 
see in the case of magnetic dots on a SC film. In that 
case, the further increase in m0 o /eo and the dot's size 
makes other vortex states more energetically favorable 
[10]. Here, the domain's size is also controlled by domain 
wall energy. Therefore, the larger the idw, the larger the 
domain size. As a result, new chain states might be more 
energetically favorable. In TV = 3 state, how the equilib- 
rium domain size, vortex-vortex distance on the chain, 
and vortex chain-domain wall distance change accord- 
ing to different values of m</>o/£o and idw are shown in 
Figs. 6, 7, 8. As seen in these figures L/A,b/A and a/A fol- 
low the same pattern as that of TV = 1 and N — 2 states. 
As number of chains per domain increases, L/A,b/A and 
a I A increase. 



5» 



- 1 



E dw = 



10 



In addition to equilibrium energies, equilibrium unit 
cell areas and vortex densities of unit cells for N = 1, 
N = 2 and N — 3 states are given in Table. IV. At fixed 
mfo/eo, unit cell area expands with further increase of 
idw However, it shrinks when m0o/£o increases. On the 
other hand, unit cell area expands with more more vor- 
tex chains at fixed values of m</>o/£o and idw except, unit 
cell areal of N = 1 state is larger than those of N = 2 
and N = 3 states, when m^o/eo = 5 and id w — 10. We 
think that N — 1 state in this case is energetically very 
close to instable region in which there is monodomain 
only. Furthermore, from our results, we expect that the 
most energetically favorable state has the highest vortex 
density. 
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FIG. 6. L/A versus m0o/£o for N = 3 state. 
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FIG. 8. a/A versus m^o/eo for N = 3 state. 
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FIG. 7. b/A versus mcpo/eo for AT = 3 state. 
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V. CONCLUSIONS 

In this article, we reported our results on the lat- 
tices of discrete vortices in stripe domains in FSB based 
on a method based on Maxwell-London equations. If 
£<iw < 4t7j 2 , the continuum approximation becomes in- 
valid. Instead, we considered the discrete lattice of vor- 
tices in which the vortices were considered to be situated 
on chains directed along the stripes. We first analyzed 
the vortex configurations up to two vortex chains. De- 
pending on the magnetization and the magnetic domain 
wall energy, the equilibrium energy, the positions of the 
vortices and the equilbrium domain size are calculated. 
According to our calculations, in equilibrium, vortices 
on the either side of the magnetic domain walls are not 
side by side on the chains; instead, they are shifted by a 
half period along the stripe, while they are side by side 
in the same domain. We also checked whether different 
vortex chain states can be energetically favorable. To do 
so, we calculated the equilibrium energy for N = 3 state 
whose equilibrium configuration was assumed to be the 
superposition of those for N = 1 and N — 2 states. The 
comporison of energies implies that more vortex chains 
can emerge spontanously in the larger domains. 

In numerical calculations, we also found that the vor- 
tex lattice is stable for Ed w > 4m 2 . At this point, the 
domain size is noticably larger than the effective pene- 
tration depth A, so the continuum approximation is valid. 
Therefore, we expect that the domain nucleation starts 
in the continuum regime. This problem is left for the fu- 
ture research. At constant £dw, with increasing m4>o/e v , 
the equilibrium size of the domain decreases. In addition, 
the vortices on the chain get closer to each other. These 
results agree with those obtained in the continuum ap- 
proximation. As £<2iu/4m 2 increases, we expect that new 
vortex chains develop within the domains. We leave the 
detailed analysis of this problem to other publication. 
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In this appendix, the detailed analysis of series is given. 
First, the series in the energy calculations of the periodic 
systems are analyzed; second, the detailed calculation of 
the vortex density is shown. The series we encounter in 
the energy calculations fall into two categories. In the 
first category, we sum over one variable. The series in 
this category are in the form of Y^r=T l/ r - Employing 
the Euler-Maclaurin summation formula [11], the sum- 
mation is found with logarithmic accuracy as 



E 



i 



In rv, 



+ C. 



(Al) 



where C <~ 0.577 is the Euler-Mascheroni constant. If 
the summation is performed over only odd integers, we 
can still transform our series to (Al). Namely, 



E 

r=0 



1 



2r + 1 



/2 



E --- E 1 



r=l 



r=l 



ln(2r m , 



+ l) + C-ln(=) 
, + C + 21n2). 



C 
2' 



(A2) 

(A3) 
(A4) 



The other double series of interest here are in the form 



of 



r— OO S—OO 



J (' T ) E E x 2 r 2 +s 2' 



(A5) 



r— — oo S— — QO 



where x is an arbitrary constant. Although (A5) is log- 
arithmically divergent, the sum over one of the vari- 
ables can be done easily. To this end, we perform 
the sum over s first. In doing so, Eq. (A5) becomes 
(2tt/x) J2T=i coth(irxr)/r [12]. This series is logarithmi- 
cally divergent. In order to get the logarithmic term , we 
add and subtract l/r. Using the result in (Al), finally 
we get 



I(x) « - 



coth(7rxr) — 1 



r=l 



+ ln 

' ma r + C 



(A6) 



Employing the same techniques, we give the results of 
the different versions of Eq. (A5) below: 



r— oo s— oo 

E ^ 



-OO S — — oo 



1 _ 2tt 

x 2 (2r + l) 2 + s 2 ~ ~~x 



E 

r=0 



coth((2r + 1)ttx) — 1 lnr. 



2r+ 1 



+ 



+ 



C 



(A7) 



r— co s— oo 

E E 



-OO S — — OQ 



1 



:r 2 (2r + l) 2 + (2s + l) 2 



g tanh((2r + l)fQ - 1 + lnr r) 



.r=0 



2r + l 



C 
+ 2 



(A8) 



In (A7) and (A8), we use J2T=o Vfe 2 + ( 2s + l ) 2 ) = wa Y, using sin 2 ((2r + l)y) = (1 - cos(2(2r + l)y))/2 or 
7rtanh(7r?//2)/(4y). In the presence of sin ((2r + l)y) or cos 2 ((2r + l)y) = (1 + cos(2(2r + l)y))/2. For example, 
cos 2 ((2r + l)y), the series can be calculated in a similar 



r— oo s— oo 

E E 

r— — oo S — — QO 



sin 2 ((2r + l)y) _ 2n 
(x 2 (2r + l) 2 + .s 2 ) ~~ ~x 

In r 



E 

r=0 



sin 2 ((2r + l)y)(coth((2r + l)?ra;) - 1) 
2r + 1 



+ 



ln|cot(y/2)| C 
4 4 



(A9) 



r— oo s— oo 

E E 

r— — oo S— — QO 



cos 2 ((2r + I)?/) _ 2tt 
(x 2 (2r + l) 2 +s 2 ) ~ ~ 



^ sin 2 ((2r + l)y)(coth((2r + l)nx) - 1) 

r=0 



2r + 1 



lnr maa; ln|cot(j//2)| C 
4 4 4 



(A10) 



We use 



E 

r=0 



cos((2r + 1)6>) _ In | cot (0/2) | 
2r + 1 ~~ 2 ' 



(All) 
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